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ADAPTIVE  MODULATION  SCHEMES  FOR  OFDM  AND 
SOQPSK  USING  ERROR  VECTOR  MAGNITUDE  (E  VM)  AND 

GODARD  DISPERSION 


Jieying  Han,  Brett  T.  Walkenhorst,  Enkuang  D.  Wang 
Georgia  Tech  Research  Institute 
Atlanta,  Georgia,  USA 


ABSTRACT 

In  this  paper,  we  develop  a  new  approach  which  enables  adaptation  across  two  modulation  schemes 
in  the  iNET  standard:  orthogonal  frequency-division  multiplexing  (OFDM)  and  shaped-offset 
quadrature  phased-shift  keying  (SOQPSK).  We  present  the  error  vector  magnitude  (EVM)  for 
OFDM  and  second-order  Godard  dispersion  (D(2))  for  SOQPSK  as  our  link  metrics  that  mea¬ 
sure  the  degradation  due  to  thennal  noise  and  channel  effects  and  then  derive  the  mathematical 
relationship  between  these  two  metrics.  This  relationship  enables  us  to  utilize  a  set  of  empirically- 
derived  rules  that  incorporate  both  modulation  schemes. 


I.  INTRODUCTION 

In  current  telemetry  systems,  a  constant  phase  modulation  (CPM)  scheme  such  as  CPM-FM  or 
SOQPSK  is  typically  used  to  transmit  data  from  an  airborne  test  article  (TA)  to  a  receiving  ground 
station  (GS).  The  chosen  scheme  will  have  a  fixed  data  rate  and  coding  rate  and  be  unable  to  adapt 
its  transmission  to  changing  channel  conditions.  In  an  attempt  to  improve  the  spectral  efficiency 
of  such  telemetry  systems,  the  Georgia  Tech  Research  Institute  (GTRI)  is  developing  an  adaptive 
modulation  and  coding  capability  called  link-dependent  adaptive  radio  (LDAR).  The  main  purpose 
of  the  algorithm  is  to  maximize  the  throughput  of  the  telemetry  link  while  maintaining  a  minimum 
level  of  reliability,  measured  as  bit  error  rate  (BER).  In  order  to  accomplish  this  objective,  the 
channel  must  be  estimated  and  link  quality  assessed  to  determine  what  mode  of  operation  would 
meet  the  objective.  Multiple  parameters  may  be  changed  in  attempting  to  maximize  throughput 
while  maintaining  link  quality,  such  as  modulation  types,  forward  error  correction  schemes  and 
rate,  signal  bandwidth,  and  number  of  sub-carriers. 

The  idea  of  adapting  modulation  and  other  parameters  to  meet  the  demands  of  the  current  channel 
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is  not  new;  there  have  been  numerous  papers  published  on  the  idea  such  as  [1,2],  and  others.  How¬ 
ever,  these  implementations  are  somewhat  limited  in  scope  compared  to  what  we  proposed  in  the 
LDAR  project,  and  very  few  studies  address  the  adaptation  of  so  many  parameters  simultaneously. 

This  paper  focuses  on  enabling  adaptation  across  two  modulation  schemes  found  in  the  iNET 
standard:  OFDM  and  SOQPSK.  Seeking  to  find  metrics  that  measure  degradation  due  to  thennal 
noise  and  channel  effects,  we  chose  error  vector  magnitude  (EVM)  as  our  metric  for  OFDM.  This 
metric  measures  the  difference  between  the  received  symbols  and  their  expected  positions  in  the 
I/Q  plane.  Because  of  the  partial  response  filter  of  SOQPSK  and  its  associated  method  of  demod¬ 
ulation  using  the  Viterbi  decoder  [3]  or  similar  decoding  method  such  as  the  BCJR’s  sequential 
decoding  scheme  [4],  EVM  has  no  meaning  for  SOQPSK  since  there  is  no  soft  decision  in  the 

demodulation,  of  the  signal.  We  therefore  chose  a  metric  called  second-order  Godard  dispersion, 
defined  as  Du>  =  E[(  \y(n)  \  )z  —  R2Y]  where  y(n)  is  the  received  signal  and  R2  is  the  amplitude 

of  the  constant  modulus  signal,  to  measure  the  quality  of  a  received  SOQPSK  signal.  This  metric 
measures  the  width  modulus  error  of  a  constant  modulus  signal  and  therefore  captures  the  effects 
of  both  thermal  noise  and  channel  effects.  In  order  to  dynamically  adapt  the  modulation  type  be¬ 
tween  these  two  schemes,  we  derive  in  this  paper  the  mathematical  relationship  between  EVM  and 
Godard  dispersion,  which  enables  us  to  utilize  a  set  of  empirically-derived  adaptation  rules  that 
incorporate  both  modulation  schemes. 

The  paper  is  organized  as  follows.  Section  II  introduces  the  SOQPSK  and  OFDM  modulation 
schemes  that  are  used  in  aeronautical  telemetry  applications.  In  Section  III,  we  present  the  link 
metrics  of  EVM  and  Godard  dispersion  (D(2)).  Section  IV  develops  the  mathematical  relationship 
between  EVM  and  Di2\  followed  by  numerical  results  in  Section  V.  Finally,  we  conclude  the  paper 
in  Section  VI  and  discuss  our  future  work. 


II,  MODULATION  SCHEMES 

Description  of  SOQPSK  The  SOQPSK  modulation  scheme  [6-8]  is  defined  as  a  continuous 
phase  modulation  (CPM)  with  baseband  modulated  signal  expressed  as 


x(t)  =  £^exp[/(<p(t  a)  +  cpo)] 

lb 


(1) 


where  Eb  is  energy  per  bit,  Tb  is  the  bit  duration,  (po  is  an  arbitrary  phase  that  will  be  set  to  0  for 
this  work,  and  the  information  carrying  phase  is  given  by 
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where  g(t)  is  the  frequency  pulse,  h  =  1/2  is  the  modulation  index,  and  an  e  {—  1 , 0, 1 }  are  the 
ternary  input  symbols,  which  are  related  to  the  binary  data  bits  bn  6  {0,  \}  by 
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Different  versions  of  SOQPSK  differ  by  their  respective  frequency  pulse.  In  this  paper,  we  use 
the  telemetry  version  SOQPSK-TG,  which  is  partial-response  with  pulse  duration  of  L  =  8  and 
a  frequency  pulse  given  by  the  product  of  a  raised-cosine  frequency-domain  window  and  time- 
domain  window  given  by  the  following  expression 


9(t)  =  C 
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According  to  the  iNET  standard  [5],  the  parameters  are  fii  =  0.7,  =  1.25,  T\ 

T2  =  0.5.  The  constant  C  is  chosen  to  make  g(t)  =  1/2  for  t  >  2(T\  +  T2)Tb- 
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1.5,  and 


The  discrete-time  SOQPSK  system  model  is  presented  in  Fig.  1,  where  xn  is  a  discrete-time  SO¬ 
QPSK  waveform,  hi  is  the  multipath  channel  impulse  response,  L  is  the  channel  length  and  wn  is 
an  additive  white  Gaussian  noise.  The  output  of  the  channel  yn  can  be  written  as 


yn  =  h  *  xn  +  wn  = 


V 

xn-ih,  +  wn. 


(6) 
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Figure  1:  SOQPSK  system  model 


Description  of  OFDM  The  OFDM  [9]  implementation  is  illustrated  in  Fig.2. 
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Figure  2:  OFDM  system  model 


At  each  OFDM  block  index  /,  a  set  of  N  QAM  symbols  =  {S('I^}N  1  is  transmitted  over  N 

k  k=  0 
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sub-carriers  and  these  frequency  components  are  converted  into  time  samples  by  perfonning  an 
inverse  DFT  (IFFT).  After  parallel  to  serial  conversion,  a  cyclic  redundancy  of  length  v  is  added  as 
a  piefix  in  such  a  way  that  x  ^  _  x^_m>  m  =  1;  2,  ...v  and  the  signal  is  then  serially  transmitted 
over  a  noisy  multipath  channel  with  channel  impulse  response  (CIR)  { hi}L~ * .  If  the  Cyclic  Prefix 
length  is  at  least  as  long  as  the  CIR  length  ( v  >  L)  [10],  the  received  signal  y(,)  can  be  expressed 
in  terms  of  at(,)  and  noise  vector  w  as 
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The  received  symbols  y(,),  ■  ■  ■  are  discarded  since  they  may  be  affected  by  ISI  from  the 

-1  - V 

prior  data  block,  and  they  are  not,  needed  to  recover  the  input  signal.  The  first  v  symbols  of  x(,) 
correspond  to  the  cyclic  prefix:  x(l  _  ('>  (/)  S')  TT7_ ,, _ _ _ 

-i  ~  x 


(/) 

N—l 


N—l; 


(8) 


where  II 


x  —  xM  .  We  therefore  write  (7)  as 

i A-v  N—v  v  ’ 

y  =  Hcircx{l)  +  tv 

is  an  N  x  N  circulant  matrix  with  an  eigenvalue  decomposition  of 

Hcirc  =  QhAQ  (9) 

where  Q  E  CNxN  is  the  normalized  Discrete  Fourier  Transform  (DFT)  matrix  (note  that  Q  is 
unitary  i.e.  QQH=  In),  and  the  diagonal  matrix  A  contains  the  eigenvalues  of  the  circulant 

(  1 

A  =  diag  Q 


matrix.  These  values  can  be  written  as 


h  l\ 
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(10) 


The  matrix  H  is  defined  to  be  a  diagonal  matrix  containing  the  channel  frequency  response  coef¬ 
ficients  along  its  main  diagonal. 


Taking  the  FFT  of  (8),  we  find 

Y{i)  =  Qy{i) 

=  Q(HcircX{l)  +  IV) 

=  Q(HcircQHS{i]+iv) 

=  Q(QhHQQhSU)+iv) 

=  HSU)  +  xv 

where  Iv=  Qiv.  After  equalization  [11],  the  estimated  symbol  S*  *  is  given  as 


S  =  H~lY(i)  =  H~l(HSu)  +  u)  =  S{l)  +  Hxiv 


(11) 


(12) 
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III.  LINK  METRICS 


Based  on  the  system  models  introduced  from  the  previous  section,  we  want  to  find  a  common  met¬ 
ric  that  will  allow  us  to  adapt  our  transmissions  over  both  of  these  modulation  schemes.  However, 
there  is  no  such  a  common  metric  that  exists  for  both  modulation  schemes,  we  therefore  chose 
EVM  as  our  metric  for  OFDM  and  Godard  dispersion  for  SOQPSK. 

In  order  to  adapt  over  both  OFDM  and  SOQPSK  with  two  distinct  metrics,  we  need  to  find  the 
relationship  between  EVM  and  Godard  dispersion.  The  two  metrics  are  detailed  in  the  following. 

Error  Vector  Magnitude  (EVM)  EVM  measures  the  deviation  of  the  received  symbols  from 
their  original  transmitted  positions  in  the  I/Q  plane.  Fig.3  shows  the  normalized  constellation  di¬ 
agram  for  QPSK  with  one  received  symbol.  It  has  been  shown  that  EVM  normalization  enables 
direct  comparison  for  a  given  average  power  level  per  symbol  between  modulation  types  (i.e., 
BPSK,  QPSK,16QAM,  64QAM,  etc)  [12]. 


Figure  3:  Normalized  constellation  diagram  for  QPSK 


EVM  is  defined  as  the  root-mean-square  (RMS)  value  of  the  difference  between  a  collection  of  re¬ 
ceived  symbols  and  the  original  transmitted  symbols.  These  differences  are  averaged  over  a  given, 
typically  large,  number  of  symbols  and  are  often  shown  as  a  percentage  of  the  average  power  per 
symbol  of  the  constellation,  or  mathematically: 


Jr,i 


dp  n  0 

EVMrmS  =  L.  _ 0  >  (13) 

C  -1  ^  9  2  a 

N  i,°’i 

s  i=  1 

where  Srj  is  the  received  symbol,  S0j  is  the  original  transmitted  symbol,  and  Ns  is  the  number  of 
symbols  over  which  EVM  is  averaged.  In  this  paper,  we  use 


EVM  =  (EVMrms)2- 


(14) 


Godard  Dispersion  Function  (Cost  Function)  Godard  dispersion  function  [13]  was  first  pro¬ 
posed  by  Godard  [14],  and  is  defined  as 

Dlp)  =  E[(|yn|p  -  Rp)2}  (15) 
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(16) 


E[W2p] 

Rp *  E[U„I2] 

where  xn  is  the  input  symbol,  yn  is  the  received  symbol,  Rp  is  a  constant  depending  only  on  the 
input  data  constellation,  p  is  an  integer.  For  p  =  2,  the  special  Godard  algorithm  was  developed  as 
the  Constant  Modulus  algorithm  (CMA)  independently  by  Treichler  and  co-workers  [15],  and  the 
CM  A  is  widely  used  for  blind  equalization.  In  this  paper,  we  use  p  =  2. 


IV,  RELATIONSHIP  BETWEEN  E VM  AND  GODARD  DISPERSION 

Based  on  the  description  of  EVM  and  Godard  dispersion  from  the  previous  section,  we  now  seek 
to  relate  the  two  metrics. 


EVM  for  OFDM  Assume  the  transmit  symbol  S(,)  is  one  OFDM  block  (with  N  sub-carriers) 

nonnalized  to  unit  energy  and  S<'l'>  denotes  the  estimated  symbol  defined  in  (12),  the  EVM  is 
defined  as  follows: 
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where  Hjj  is  the  jth  diagonal  element  of  H  (defined  in  (10)) .  As  iv=  Qiv,  we  have 
E[tixtf]  =  E[(Qiv)(Qiv)h]  =  E[QivwhQh]  =  Qo2IQH  =  o2I 

W  W 

Thus  u)  ~  CN  (0,  o2 1 )  and  o2  =  <j2 .  Therefore, 
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Godard  Dispersion  for  SOQPSK  Let  xn  denote  the  transmitted  complex  signal,  yn  denote  the 
received  signal,  h  denote  the  channel  impulse  response  of  length  L,  and  wn  denote  the  noise  with 
distribution  w  ™  CN  (0,  o^I).  The  received  signal  is  then  given  by  yn  =  h  *  xn  +  wn.  Assume 
we  normalize  the  transmitted  signal  to  unit  energy:  E[  \xn  | 2]  =  (72  =  1 ,  and  the  channel  has  unit 
energy:  q  \  hi\ 2  =  1  •  The  second-order  Godard  dispersion  (D<2))  for  SOQPSK  is  defined  as: 


l(  \2I 

D(2)  =  E  | yn 1 2  -  1 


(20) 


We  model  xn  as  a  complex  variable  with  zero  mean  and  variance  o2 .  Given  E[(xn)2]  =  E[(x*)2]  = 

x  n 

E[xn\  =  0,  we  find  the  following  (see  Appendix  for  the  complete  derivation): 

D(2)  =  2  \hi\2\hp\2-  |/7,|4  +  2ct4+2o2-1  (21) 

w  w 
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Mapping  Between  EVM  and  Godard  Dispersion  From  the  previous  sections,  we  have  ex¬ 
pressed  EVM  and  second-order  Godard  Dispersion  in  tenns  of  the  noise  variance  and  channel 
impulse  response.  We  now  relate  these  metrics  to  one  another  by  solving  for  o2w  in  (19)  and  sub¬ 
stituting  into  (21)  to  obtain 


\\i 


k:l 


D{2)  =  2  \hi\2\hp\2  —  |  h,\‘ 

1=0  p=0  1=0 


I  ( _  EVM  V 


1  f  >  N 

N  jj= 


EVM  | 

-l 


1  (22) 


ii«SV 


1  N 
N  jj=  1 


I "JJ  I2 


where  N  is  the  number  of  sub-carriers,  /?,  is  channel  impulse  response  and  Hjj  is  the  frequency- 
domain  channel  response  of  the  fh  subcarrier. 


V.  SIMULATION  RESULTS 

We  now  perform  simulations  to  compare  our  theoretical  mapping  to  metrics  computed  from  time- 
domain  simulated  data  using  an  AWGN  channel.  The  theoretical  data  is  calculated  using  (19)  and 
(21),  the  simulated  data  is  calculated  using  (13)  and  (15)  with  simulated  SOQPSK  and  OFDM 
waveforms.  The  simulation  results  were  obtained  using  2/3  low  density  parity  check  (LDPC)  code 
rate.  The  EVM  values  were  calculated  using  OFDM/QPSK  with  64  sub-carriers.  Both  EVM  and 
D(2)  were  simulated  using  100  iNET  bursts,  each  burst  consists  of  49152  coded  bits  (8  LDPC 
codeblocks). 

Table  1  presents  our  results  using  SNR  =  —5  dB,  0  dB  and  15  dB,  the  simulated  data  of  EVM 
and  D(2)  are  both  very  close  to  the  theoretical  data,  which  verifies  our  mathematical  expressions  in 
(19)  and  (21).  The  “mapped  D(2)  from  EVM”  is  essentially  identical  to  the  theoretical  D(2*,  which 
verifies  our  mathematical  mapping  in  (22). 

More  simulation  results  are  shown  in  Fig.4  with  SNR  ranges  from  —5  dB  to  20  dB.  The  top  figure 

is  for  EVM  while  the  bottom  is  for  Godard  dispersion  (D(2)).  It  can  be  seen  from  the  figure  that  the 
simulated  EVM  and  D(2)  are  both  very  close  to  the  theoretical  data,  the  average  percentage  error 
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between  “mapped  D(2)  from  EVM”  and  “theoretical  D(2)”  is  only  about  0.0682%. 
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Table  1:  Simulation  results  comparing  to  theoretical  results 


SNR(dB) 

EVM 

Godard  Dispersion  (D(Z}) 

Theoretical 

Simulated 

Theoretical 

Simulated 

Mapped  from  EVM 

-5 

3.1623 

3.1642 

26.3246 

26.3318 

26.3525 

0 

1.0000 

0.9993 

4.0000 

3.9968 

3.9960 

15 

0.0316 

0.0316 

0.0652 

0.0652 

0.0652 

Figure  4:  Simulation  results  comparing  to  theoretical  results 


VI.  CONCLUSION 

We  presented  a  theoretical  mapping  between  two  metrics  to  enable  the  adaptation  across  two  mod¬ 
ulation  schemes  in  the  iNET  standard:  OFDM  and  SOQPSK.  We  derived  the  mathematical  map¬ 
ping  from  EVM  to  second-order  Godard  dispersion,  which  enables  us  to  utilize  a  set  of  empirical 
adaptation  rules  that  incorporate  both  modulation  schemes  [16].  Moreover,  we  verified  our  map¬ 
ping  using  an  AWGN  channel  with  experimental  results  very  close  to  theoretical  results.  Future 
work  will  integrate  this  mapping  to  our  adaptation  rules  in  [16]  and  test  using  representative  chan¬ 
nel  models  we  have  developed;  then  we  will  implement  these  rules  in  hardware  and  test  in  real 
wireless  channels. 
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Appendix:  Derivation  of  (21) 


(  \2  (  \2 
|yn|  -1  =  \h*xn  +  w„\z-'\ 


( L\1  \  V 

=  (  Xn-ih,  +  Wn  )(  Xn-ihi  +  Wny  -  1 

/=o  /=0  _  _ 

L-1L-1  ZSl  lNL 

VXX  ih  x*  h*  +  X*  +  X„/^W*  +  |Wnl 

=  Xn—  I  n—k  k  n—l  /  —  In 

1=0  k=  0  /=0  /=0 

/  n-k^ k^n—ph pXn_qhq 

1=0  k=0  p=  0  q=0 

XX 


\ 

1  2 


/=0  k=  0 
L—lL—1 


1  n-A  WnXn_khkWn 


w 

P  n 


X  X  ihw*x  h  w* 

+  xn_  1  n  n~k  k  n 

1=0  k=0 

+  (|w„|2-1)2 

L-l  L-l  L-l 

XXX  ih  x*  u*x*  h*wn 

+  2  Xn_  1  n-k  k  n-p  p 

1=0  k=  0  p=0 

ihx*  h*xn  ph 

+  2  Xn“  1  n-k  k 

1=0  k=  0  p=  0 

+  2(  |  wn|2  —  1 )  xn  ih  X*  h* 

—  I  n-k  k 

1=0  k=0 

XlXi 

+  2  Xn-lh>nXn-khktfn 

1=0  k=  0 

-i-2(|  wn  | 2  —  1 )  x*n_i  hjwn 

+  2( | wn|2  —  1 )  xn  ,hw* 

—  I  n 

1=0 


Using  the  equation  labels,  we  have 

l( 

D(2)=E  |  yn  1 2  -  1 


\2' 


(23) 

(24) 

(25) 

(26) 

(27) 

(28) 

(29) 

(30) 

(31) 

(32) 


=  E[(23)]  +  E[(24)]  +  E[(25)]  +  E[(26)]  +  E[(27)]  +  E[(28)] 
+  E[(29)]  +  E[(30)]  +  E[(31 )]  +  E[(32)] 
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Approaching  these  in  turn,  we  have 


E[(  23)]  = 


1=0  p=0 

( l=k,p=q,l/=p 


L- 1 L- 1 


|hp|2E  |xn  _/|  2  E  | Xn-p 


I  II  I 

|/7/|2|^|2E  |xn_,|  2  E  \xn-k\2 


(l=q,k=p,l/=k 


L_1  I  I 

^  |h,|4E|xn_/|4  +0 

+  ,  „  (oi 


E(  24) 

I 

E  (25) 

I 

E  (26) 

e'(27) 

I 

E  (28) 

I 

E  (29) 

e'(30) 

I 

£(31) 


+  (others) 

( l=k=p=q ) 

W  K1 

=  2  |ft,l2IM2-  |ft,l4 

/=0  p=0  /=0 

\  I  II  II  I 

=  E  OC,)2  E  (A?;)2  E  {Wnf  =  0 

1=0 

I  II  II  I 

=  E  (xn_,)2  E  {hi)2  E  «)2  =0 

=f  I 

=  E  |  wn  | 4  +  1  —  2 1 | 2  =2cr4+1-2cr2 

1/1/  1/1/ 

L-l  L-t  L-l  . 

I  *  x*  E  h/h*h*  E  wn  =  0 

=  2  E  Xn~  n-k  n-p  k  p 

1=0  k=0  p=0 

1Sl1SJ> i.  I  II  III 

=  2  E  xn  *  Xn  p  E  h  h*hp  E  w*  =  0 

—  n-k  -  /  k  n 


1=0  k=o  p= 0 

^  '  2 1 

=  2(a2-1)  \h,\2E  \xn-i\  =  2{o2w—  1 ) 

/=°  . 

rV/N.  I  TV 

=  2  E\  Xn  k\  \h*hk\  I wnw*  =  2  a2a2 \h,\2  =  2a2 

yfi—l  -[=11=  n  X  w  w 

1=0  k=0  1=0 

I  k-i  |  III 


1=0  k=0 

I 


-2 E  i/i/n(|i/i/n|2-1) 


*  * 

Eh,  =0 

/=0  Z- 

I  LI  I  III 


E  (32)  -2E  i/i/*(|i/i/n|2- 1)  ^  xn-i  E  hi  -  0 

/=o  £ 

Finally,  we  have  D(2)  as  follows: 


D(2)  =  2^^  | h, | 2 1 hp | 2  - h, \ 4  +  2cr4  +  2CJ2  -  1 

w  w 

1=0  p=  0  1=0 


11 


